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We present three different matrix bases that can be used to decompose density 
matrices of d-dimensional quantum systems, so-called qudits: the generalized Gell- 
Mann matrix basis, the polarization operator basis, and the Weyl operator basis. 
Such a decomposition can be identified with a vector — the Bloch vector, i.e. a 
generalization of the well known qubit case — and is a convenient expression for 
comparison with measurable quantities and for explicit calculations avoiding the 
handling of large matrices. We consider the important case of an isotropic two- 
qudit state and decompose it according to each basis. Investigating the geometry of 
entanglement of special parameterized two-qubit and two-qutrit states, in particular 
we calculate the Hilbert-Schmidt measure of entanglement, we find that the Weyl 
operator basis is the optimal choice since it is closely connected to the entanglement 
of the considered states. 



I. INTRODUCTION 

The state of a ci-dimensional quantum system — a qudit — is usually described by a 
d X d density matrix. For high dimensions, where the matrices become large (for composite 
systems of n particles the matrices are of even much larger dimension x rf"), a simple 
way to express density matrices is of great interest. 

Since the space of matrices is a vector space, there exist bases of matrices which can be 
used to decompose any matrix. For qubits such a basis contains the three Pauli matrices, 
accordingly, a density matrix can be expressed by a 3-dimensional vector, the Bloch vector 
and any such vector has to lie within the so-called Bloch ball [H, 0|- Unique for qubits is the 
fact that any point on the sphere, Bloch sphere, and inside the ball corresponds to a physical 
state, i.e. a density matrix. The pure states lie on the sphere and the mixed ones inside. 

In higher dimensions there exist different matrix bases that can be used to express qudits 
as (d^ — l)-dimensional vectors as well. Different to the qubit case, however, is that the 
map induced is not bijective: not every point on the "Bloch sphere" in dimensions d^ — 1 
corresponds to a physical state. Nevertheless the vectors are often also called "Bloch vectors" 
(see in this context, e.g., Refs. 0, 0, H, 0, 01)- 

In this paper we want to present and compare three different matrix bases for a Bloch 
vector decomposition of qudits. In Sec. [Ill we propose the properties of any matrix basis for 
using it as a "practical" decomposition of density matrices and recall the general notation 
of Bloch vectors. In Sees. IIIII - |V] we offer three different matrix bases: the generalized 
Gell-Mann matrix basis, the polarization operator basis, and the Weyl operator basis. For all 
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these bases we give examples in the dimensions of our interest and present the different Bloch 
vector decompositions of an arbitrary density matrix in the standard matrix notation. Next 
in Sec. |VT1 by constructing tensor products of states we study the isotropic two-qudit state 
and present the results for the three matrix decompositions, i.e. for the three different Bloch 
vectors. In Sec. IVIII we focus on the geometry of entanglement of the isotropic two-qudit 
state and calculate the Hilbert-Schmidt measure of entanglement (see, e.g., Refs. @MM, 



HI]). Its connection to the optimal entanglement witness is shown, which is determined in 



terms of the three matrix bases. Furthermore, we calculate explicitly the Hilbert-Schmidt 
measure for so-called two-parameter states, which define a plane in the Hilbert-Schmidt 
space and provide a picture for the geometry of entanglement of states in higher dimensions, 
in particular we studied the two-qubit (4-dimensions) and two-qutrit states (9-dimensions). 
The mathematical and physical advantages/disadvantages by using the three different matrix 
bases are discussed in Sec. IVIIIt where also the final conclusions are drawn. 



II. PRELIMINARIES 



A qudit state is represented by a density operator in the Hilbert-Schmidt space acting 
on the d-dimensional Hilbert space T-C^ that can be written as a matrix — the density 
matrix — in the standard basis {\k)} , with k = 1,2, . . . d or k = 0,1,2, . . . d — 1. 

Properties of a "practical" matrix basis. For practical reasons the general properties of a 
matrix basis which is used for the Bloch vector decomposition of qudits are the following: 

i) The basis includes the identity matrix 1 and d — 1 matrices {Ai} of dimension d x d 
which are traceless, i.e. TrAi = . 

ii) The matrices of any basis {Ai} are orthogonal, i.e. 

TiAlAj = NSij with NeR. (1) 



Bloch vector expansion of a density matrix. Since any matrix in the Hilbert-Schmidt space 
of dimension d can be decomposed with a matrix basis {Ai}, we can of course decompose a 
qudit density matrix as well and get the Bloch vector expansion of the density matrix, 

p = + (2) 

where 6 ■ F is a linear combination of all matrices {Ai} and the vector b G R''^^^ with 
hi = (Fj) = TrpFj is called Bloch vector. The term ^1 is fixed because of condition Trp = 1. 



Remark. Note that a given density matrix p can always be decomposed into a Bloch 
vector, but not any vector a that is of the form ([2]) is automatically a density matrix, even 
if it satisfies the conditions Tr cr = 1 and Tro"^ < 1 since generally it does not imply a > 0. 

Each different matrix basis induces a different Bloch vector lying within a Bloch hyper- 
sphere where, however, not every point of the hypersphere corresponds to a physical state 
(with p > 0); these points are excluded (holes). The geometric character of the Bloch space 
in higher dimensions turns out to be quite complicated and is still of great interest (see 

Refs. BHBBH). 
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All different Bloch hyperballs are isomorphic since they correspond to the same den- 
sity matrix p. The interesting question is which Bloch hyperball — which matrix basis — is 
optimal for a specific purpose, like the calculation of the entanglement degree or the deter- 
mination of the geometry of the Hilbert space or the comparison with measurable quantities. 



III. THE GENERALIZED GELL-MANN MATRIX BASIS 

A. Definition and example 

The generalized Gell-Mann matrices (GGM) are higher-dimensional extensions of the 
Pauh matrices (for qubits) and the Gell-Mann matrices (for qutrits), they are the standard 
SU(N) generators (in our case N = d). They are defined as three different types of matrices 
and for simplicity we use here the operator notation; then the density matrices follow by 
simply writing the operators in the standard basis (see, e.g. Refs. 0,113): 

i) ^i^tli symmetric GGM 

Af = \3){k\ + |A;)(j|, l<3<k<d, (3) 



d(d-l) 



antisymmetric GGM 



i\3){k\+i\k){3\, l<j<k<d, (4) 



+ + , l<l<d-l. (5) 



In total we have d"^ — 1 GGM; it follows from the definitions that all GGM are Hermitian 
and traceless. They are orthogonal and form a basis, the generalized Gell-Mann matrix 
Basis (GGB). A proof for the orthogonality of GGB we present in the Appendix I A 1[ 

Examples. Let us recall the case of dimension 3, the 8 Gell-Mann matrices (for a repre- 




sentation see, e.g., Refs. [111. Il3l|) 



i) 3 symmetric Gell-Mann matrices 

= ( ? 1 , Ai^ = ( 1 , Xf 






(6) 



ii) 3 antisymmetric Gell-Mann matrices 
/ -i 0\ 

= ^ , Ai3 = ' n n n 1 \23 










—i 











i 










(7) 
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ii) 2 diagonal Gell-Mann matrices 





(8) 



To see how they generahze for higher dimensions we show the case we need for qudits of 
dimension d = 4: 
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(10) 



iii) 3 diagonal GGM 



Ai = 



/ 1 0\ 

0-100 
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A2 = ^ 
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(11) 



Using the GGB we obtain, in general, the following Bloch vector expansion of a density 
matrix 

p ^ ll + b-A, (12) 



with the Bloch vector b = {{bi''}, {b^''}, , where the components are ordered and for the 
indices we have the restrictions I < j < k < d and 1 < I < d—1 . The components are given 
by bi'' = TrA^'^p , b^'' = TrA^'^p and 6' = TrA'p . All Bloch vectors he within a hypersphere 
of radius |&| < ^/(d — l)/2d. For example, for qutrits the Bloch vector components are 
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b = (pl'^ , bl^ , bf , bl^ , , bl^ , , b'^) corresponding to the Gell-Mann matrices dS]), ([TD, ([8]) 
and |6| < . 

— * 

As already mentioned the allowed range of b is restricted. It has an interesting geo- 
metric structure which has been calculated analytically for the case of qutrits by studying 

2- dimensional planes in the 8-dimensional Bloch space or numerically by considering 

3- dimensional cross-sections In any case, pure states lie on the surface and the mixed 
ones inside. 



B. Standard matrix basis expansion by GGB 

The standard matrices are simply the d x d matrices that have only one entry 1 and the 
other entries and form an orthonormal basis of the Hilbert-Schmidt space. We write these 
matrices shortly as operators 

\j){k\, with j,k = l,...,d. (13) 

Any matrix can easily be decomposed into a "vector" via a certain linear combination of 
the matrices f|T3l) . Knowing the expansion of matrices f|T3|) into GGB we can therefore find 
the decomposition of any matrix in terms of the GGB. 

We find the following expansion of standard matrices f[T3|) into GGB : 



\j){k\ 



\ (Af + zAf ) for 3 <k 

\ {A'^J - ih^/) for ]>k 

d-j-l 

A^-i + V ^ A-'+" + i 1 for J = A; . 

2^ „=0 V^U+n)U+n+l) d 



(14) 



Proof. The first two cases can be easily verified. 

To show the last case we first set up a recurrence relation for |/)(/|, which we obtain by 
ehminating the term in the two expressions ([S]) for A' and A'~^ 



I0(/| = -\/^^ + V^"^ + + + (15) 



and we consider the case I + 1 = d 



\d-l){d-l\ = -xI4tA^^''" + xhrJ-r,^'-' + \d){d\. (16) 



2{d- 1) V 2{d- 1) 



From A'^ ^ given by Eq. ([5]) 



A'"' = i/tT^ I > :U>{jI - (i-mW\ I , (17) 



we get the Bloch vector decomposition of \d){d\ 



1/ ./(rf-lM.d-i 



\d){d\ = \ ' (18) 
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where we have apphed X]j=i = 1 — \d){d\ . 

Inserting now decomposition (fT8l) into relation (fT6l) we gain the Bloch vector expansion 
for \d — l){d — 1| and recurrence relation (1151) provides \d — 2){d — 2| and so forth. Thus 
finally we find 



\d-n)(d-n\ = -, r ^ 1a^-"-i + \" ^ \<i-n+k,l^ 

\l "^{d-n) ^ ^2{d-n + k + l){d-n + k) d' 

(19) 

the relation we had to prove, where d — n = j . □ 



IV. THE POLARIZATION OPERATOR BASIS 

A. Definition and examples 

Definition. The polarization operators in the Hilbert-Schmidt space of dimension d are 
defined as the following d x d matrices [3, Il3] : 



d 



Tlm = \r^;^Y.^sZLM\k){l\. (20) 

kA=l 



The used indices have the properties 



„ _ d-l 
^ 2 ' 



L = 0,1, ... ,2s, 
M = -L,-L + 1,...,L-1,L, 
mi = s, 1712 = s — 1, . . . , rrid = —s . (21) 

The coefficients C^™* r m are identified with the usual Clebsch-Gordan coefficients Cj"l. ^ 



of the angular momentum theory and are displayed explicitly in tables, e.g., in Ref. [14 



For L = M = the polarization operator is proportional to the identity matrix 0, 14 

Too = ^ 1 . (22) 

It is shown in Ref. [3] that all polarization operators (except Too) are traceless, in general 
not Hermitian), and that orthogonality relation is satisfied 

"^^Tl^j^^j^Ti^Mi = ^LiLa^MiMz • (23) 

Therefore the d"^ polarization operators ( l20i) form an orthonormal matrix basis — the 
polarization operator basis (FOB) — of the Hilbert-Schmidt space of dimension d. 

Examples. The simplest example is of dimension 2, the qubit. For a qubit the FOB is 
given by the following matrices {s = 1/2; L = 0,1; M = —1, 0, 1) 

^ _ 1 / 1 0\ „ fO 1 
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For the next higher dimension d = 3 (s = 1), the case of qutrits, we get 9 polarization 
operators T^m with L = 0, 1, 2 and M = —L, ...,L . According to definition (l20l) they are as 
follows (Too = 75 1) • 







^0 
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T22= 1 
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-1 








r2_i = ^ I 1 I , T2_2 = I I . (25) 

Then the decomposition of any density matrix into a Bloch vector by using the FOB has, 
in general, the following form : 

P = + J2 Yl ^lmTlm = -^l + b-f, (26) 

L=l M=-L 

with the Bloch vector h = (6i_i, 6io, 62-2, &2-1, •••5 &lm), where the components are 
ordered and given by Blm = '^'^T'ImP ■ general the components hiM are complex since the 
polarization operators Tlm are not Hermitian. All Bloch vectors lie within a hypersphere of 
radius \h\ < ^/{d — l)/d . 

In 2 dimensions the Bloch vector b = (&1-1, &107 &11) is limited by \b\ < and forms a 
spheroid [3], the pure states occupy the surface and the mixed ones lie in the volume. This 
decomposition is fully equivalent to the standard description of Bloch vectors with Fauli 
matrices. 

In higher dimensions, however, the structure of the allowed range of b (due to the posi- 
tivity requirement p > 0) is quite complicated, as can be seen already for 0? = 3 (for details 
see Ref. Nevertheless, pure states are on the surface, mixed ones lie within the volume 
and the maximal mixed one corresponds to \b\ = 0, thus |6| is a kind of measure for the 
mixedness of a quantum state. 



B. Standard matrix basis expansion by POB 



The standard matrices f|T3l) can be expanded by the FOB as follows 14 



1^) 0- 1 = E E v 1^ LM tlm . (27) 



, 2s + l 



Note that "YIim actually fixed by the condition rrtj + M 



rrii 



8 



Proof. Inserting definition (l20l) on the right-hand side (RHS) of equation (1271) we find 

RHS = ^ I ^ 2s + I ^smi,LM I 1^) (^1 = 

k,l \ L ) 

k,l 

= mj\, (28) 



where we used the sum rule for Clebsch-Gordan coefficients 14 



26+1 

c,7 



V. WEYL OPERATOR BASIS 
A. Definition and example 



Finally we want to discuss a basis of the Hilbert-Schmidt space of dimension d that 
consists of the following (f operators 



d-l 

2-ni 



Unm = ^ '"'\k){{k + m) mod d\ n, m = 0, 1, . . . , d - 1 , (30) 



k=0 



where we use the standard basis of the Hilbert space. 

These operators have been introduced in the context of quantum teleportation of qudit 
states [15| and are often called Weyl operators in the literature (see e.g. Refs. 0,[l3)- The 



d"^ operators fl30|) are unitary and form an orthonormal basis of the Hilbert-Schmidt space 
(a proof is presented in Appendix lA 31) - the Weyl operator basis fWOB). They can be used 



to create a basis of c?^ maximally entangled qudit states 16|, UM, Il9 |. 



Clearly the operator Uqq represents the identity f/oo = 1 . 

Example. Let us show the example of dimension 3, the qutrit case. There the Weyl 
operators (130|) have the following matrix form 



(31) 



U- 
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Using the WOB we can decompose quite generally any density matrix into a Bloch vector 

^ d-i 1 ^ ^ 

P = -:1 + ^nrrJJnm = 1 + b ■ U , (32) 



a ^-^ a 

n,m=0 

with n,m = 0,1, d — 1 (600 = 0). The components of the Bloch vector b = {{bnm}) 
are ordered and given by bnm = Tr Unm p ■ In general the components bnm are complex 
since the Weyl operators are not Hermitian and the complex conjugates fulfil the relation 
^nm — ^^^^b-n-m , which follows casily from definition (130|) together with the hermiticity 
of p. 

All Bloch vectors lie within a hypersphere of radius |6| < ^/ d — 1/d . For example, for 
qutrits the Bloch vector is expressed by 6 = (601, ^11, ^22 ) and |6| < V2/3. 

In 3 and higher dimensions the allowed range of the Bloch vector is quite restricted within 
the hypersphere and the detailed structure is not known yet. 

Note that in 2 dimensions the WOB as well as the GGB coincides with the Pauli matrix 
basis and the FOB represents a rotated Fauli basis (where cr± = ^ (cxi ± ia2)), in particular 

{Uoo,Uoi,Uio,Uu} = {1,0-1,0-3,^(72}, (33) 

{l,Af,Af,Al} = {1, (71,(72, (73} , (34) 

{Too,Tii,Tio,Ti_i} = |— 1, -0-+, — 0-3, o-_| . (35) 

B. Standard matrix basis expansion by WOB 

The standard matrices ( fT3l) can be expressed by the WOB in the following way 

= ■^Y.''''^''^nk-j)modd. (36) 

Proof. We insert the definition of the Weyl operators (l30i) on the right-hand side (RHS) 
of Eq. dSSD, use Eq. flA24D and get 

RHS = - ^e^'(^-^')|r)((r + fc-j)modrf| 

l,r=0 

^ d~l d-1 

= \j){k\ + Yl Xl^^'^''''^ 

r^j, r=0 1=0 

= \j){k\. □ (37) 



VI. ISOTROPIC TWO QUDIT STATE 

Now we consider bipartite systems in a. d x d dimensional Hilbert space ^Tig. The 
observables acting in the subsystems T-La and TIb are usually called Alice and Bob in quantum 
communication. 
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Quite generally, a density matrix of a two-qudit state acting on can be decom- 

posed in the following way (neglecting the reference to A and B) 

p = ^ 1 ® 1 + nj Tj (g) 1 + mj 1 (g) r, + (g) Tj , rii, rrii, Cij G C , (38) 

where {Fj} represents some basis in the subspace T-C^ . The term Cij Fj (g F^ always can be 
diagonalized by two independent orthogonal transformations on Fj and Tj [20]. Altogether 
there are (d^)^ — 1 independent terms. 

However, for isotropic two-qudit states — the case we consider in our paper — the second 
and third term in expression (!38l) vanish and the fourth term reduces to cu F j ® F j , which 
implies the vanishing of (c/^ — l)^ + ((i^ — 1) = (i^((i^ — 1) terms. Consequently, for an 
isotropic two-qudit density matrix there remain d"^ — 1 independent terms, which provides 
the dimension of the corresponding Bloch vector. Thus the isotropic two-qudit Bloch vector 
is of the same dimension — lives in the same subspace — as the one-qudit vector, which is a 
comfortable simplification. 

Explicitly, the isotropic two-qudit state p^a'' is defined as follows 

pL'^ = ak+>(0t| + aeR, -^^<a<l, (39) 

where the range of a is determined by the positivity of the state. The state | a Bell 
state, is maximally entangled and given by 

14) = 4^ E i-^') ® 1^') ' (40) 

j 

where denotes the standard basis of the d-dimensional Hilbert space. 



21 



23 



A. Expansion into GGB 

Let us first calculate the Bloch vector notation for the Bell state |(/)^|_) {(f)'^\ in the GGB. 
It is convenient to split the state into two parts 

= ^E i^')(^i®i^')(^i 
j,k=i 

= A + B, (41) 

where A and B are defined by 

A := \ j2\j){k\®\j){k\+^-J2\k){j\®\k){j\, (42) 

j<k j<k 

B := ^El^')01®b)01, (43) 



and to calculate the two terms separately. 
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For term A we use the standard matrix expansion (fT^ for the case j ^ k and get 



.j<k 



j<k 



1 J2 (Af ® Af - Af ® Af ) . 



(44) 



For term B we need the case j = A; in expansion ( |T^ and obtain after some calculations 
(the details are presented in Appendix IA2p 



d-l 



B 



l5^A™®A- + ll®l 

m=l 



(45) 



Thus all together we find the following GGB Bloch vector notations, for the Bell state fl4T|) 

= 4l®l + ;^A, (46) 



and for the isotropic two-qudit state (!39|l 

1 

rf2 



where we defined 



d-l 



A := ^ Af ® Af - Af ® Af + ^ A" ® A'" . 



(47) 



(48) 



i<j 



i<j 



m=l 



B. Expansion into POB 

Now we calculate the Bell state \4>'\.) in the POB. Using expansion ( |27l) and the sum 
rule for the Clebsch-Gordan coefficients 14 



we obtain 



9r 4- 1 



a,7 



26 + 1 



(49) 



1 



1 ^ y/{2L + l){2L' + l) ( ^ \ 

^ 2s + 1 I ^ ^smj,LM^sm.j,L'M I 

l^ v/(2L + l)(2L- + l) 



Tlm ® T^iM 



2L + 1 

L,L' 



(50) 
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where we extracted the unity (recall Eq. (l22l) ) and defined 

T ■■= Tlm® Tlm. (51) 

L,M^O,0 

Result (150|1 provides the POB Bloch vector notation of the isotropic two-qudit state (l39ll 

p(f) = jjl®l + ^T. (52) 

C. Expansion into WOB 

Finally we present the Bell state in the WOB (the details for our approach using the 
standard matrix expression (15^ can be found in the Appendix lA 41 see also Ref. [3]) 

with 

d-l 

U := Y,Uira®U.i^, (/,m)^(0,0), (54) 

/,m=0 

where negative values of the index / have to be considered as mod d , and from formula (|53|) 
we find the WOB Bloch vector notation of the isotropic two-qudit state 

= + (55) 

VII. HILBERT SCHMIDT MEASURE APPLICATIONS OF THE MATRIX 

BASES 

A. Entangled isotropic two qudit states 



In Ref. [ll| the connection between the Hilbert-Schmidt (HS) measure of entanglement 
js], [q], [i3| and the optimal entanglement witness is investigated. Explicit calculations for 
both quantities are presented in case of isotropic qutrit states. For higher dimensions, the 
isotropic two-qudit states, the above quantities are determined as well but in terms of a 
rather general matrix basis decomposition. With the results of the present paper we can 
calculate all quantities explicitly. Let us recall the basic notations we need. 

The HS measure is defined as the minimal HS distance of an entangled state Pent to the 
set of separable states S 

DiPent) ■= min||p-pent|| = ||Po-Pcnt|| , (56) 

where po denotes the nearest separable state, the minimum of the HS distance. 

An entanglement witness A E A {A = Aa<^Ab , the HS space of operators acting on the 
Hilbert space of states] is a Hermitian operator that "detects" the entanglement of a state 
Pent via inequalities [ly, 24, [2I, 26] 



(Pent,^) = Trpent^ < 0, 

{p,A) = TtpA > WpeS. (57) 
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An entanglement witness is "optimal" , denoted by Aopt , if apart from Eq. (IFH) there exists 
a separable state po E S such that 

(po,Apt) = 0. (58) 

The operator Aopt defines a tangent plane to the set of separable states S and all states Pp 
with {pp, Aopt) = lie within that plane; see Fig. [TJ 



According to Ref. [10|, we call the lower one of the inequalities fl57l) a generalized Bell 
inequality, short GBI. "Generalized" means that it detects entanglement and not just non- 
locality. Re-writing Eq. fl57|l as 



{p,A) - {Pent, A) > ypeS, (59) 
the maximal violation of the GBI is defined by 

-B(pent) = max I min (p, A) - (pent, A) I , (60) 

A, ||A-al||<l y pes y 

where the maximum is taken over all possible entanglement witnesses A, suitably normalized. 

Then an interesting connection between the HS measure and the concept of entanglement 
witnesses is given by the Bertlmann-Narnhofer-Thirring Theorem, illustrated in Fig. [1] [13] . 



Theorem 1. 

i) The maximal violation of the GBI is equal to the minimal distance of pent to the set S 

^(Pent) = 5(pent) • (61) 

a) The maximal violation of the GBI is attained for an optimal entanglement witness 

^ ^ ^ Po - Pent - (PO,PO - Pent) 1 ^g2) 
IIPO — Pent II 
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Thus the calculation of the optimal entanglement witness Aopt to a given entangled 
state Pent reduces to the determination of the nearest separable state po • In special cases 
Po might be easy to find but in general its detection is quite a difficult task. 

Now let us apply the matrix bases we discussed in the previous sections and calculate 
the quantities introduced above. As an entangled state we consider the isotropic two-qudit 
state p^a^''^^^, that is the state pL'^'* ( !39l) for < a < 1. 

Starting with the GGB we can express that state in our Bloch vector notation by formula 
(j47|l . Then the nearest separable state is reached at a = 

pS^'-pH^- J.l«l + ^^A. (63) 



It provides the HS measure 



and the optimal entanglement witness (|62 



Aop.(p«„.) = y^I»I-5;^A. (65) 



where we used the HS norm ||A|| = 2^/3? — 1 . 

Clearly, the maximal violation of the GBI B equals the HS measure D 

r>( id.) \ _ _ I id) A 
[Pa, ent) — \ Pa, ent ; ^opt 



^^"a--v^) = (66) 



d V d+l 

For expressing above quantities by the matrix bases FOB and WOB it suffices to calculate 
the proportionality factors between A, T and U . By comparison of the three forms for the 
isotropic qudit state (H7I) . (132]) and (ITO]) we find 

A = 2T and T = -U . (67) 

a 

It provides the following expressions, for the FOB 



AoAPaU) = ] \ ^ 1 ® 1 - T , (69) 



and for the WOB 



d\ d+l Vd^^ 



Po^ = Pal^ = + 3^7:^^' (70) 



Aopt(p?lj-i^^l^l-^^f/. (71) 

Of course, the HS measure D{p^^\^^) remains the same expression flMIl independent of the 
chosen matrix basis, which can easily be verified using ||T|| = — 1 and \\U\\ = d\/ d^ — \ . 
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B. Two parameter entangled states — qubits 

As an application of our Bloch vector notation in the several matrix bases we want to 
determine the HS measure of entanglement for the following two-qubit states which are a 
particular mixture of the Bell states \<P^) 

Pc, = 1 + + f + . (72) 

The states fl72]) are characterized by the two parameters a and (3 and we will refer to the 
states as the two-parameter states. Of course, the positivity requirement pa,/3 > constrains 
the possible values of a and /3, namely 

1 1 

a < -f3 + l, a > -/3--, a</3 + l, (73) 

3 3 

which geometrically corresponds to a triangle, see Fig. [2l 

According to Peres [13] and the Horodeckis [l^] the separability of the states is determined 
by the positive partial transposition criterion (PPT), at least in dimensions 2®2 and 2®3 . 
States fl72p which are positive under partial transposition have the following constraints 

a > (3-1, a < i/? + 1, a>-(3-l, (74) 

and correspond to the rotated triangle; then the overlap, a rhombus, represents the separable 
states, see Fig. [21 

In the picture drawn in Fig. [2] the orthogonal lines are indeed orthogonal in HS space. 
Therefore the coordinate axes for the parameter a and (3 are necessarily non-orthogonal. In 
particular, the a axis has to be orthogonal to the boundary line a = —(3 — 1, and the (3 axis 
has to be orthogonal to a = /3 + 1. 

The two-parameter states pa,i3 define a plane in the HS space. It is quite illustrative to see 
how this plane is located in the 3-dimensional spin space of the density matrices, where the 



4 Bell states form a tetrahedron due to the positivity condition |10l . |28| . |29| . Applying PPT 
the tetrahedron is rotated producing an intersection — a double pyramid — which represents 
the separable states. This is shown in Fig. [31 

To calculate the HS measure (13^ for the two-parameter qubit state f[721) we express the 
state in terms of the Pauli matrix basis, which is indeed just the GGB and equivalent to the 
WOB for dimension d = 2 (see Sec. IV Ap 



Pa,p = ^ (1 + a (o-i (g) (Ji - 0-2 ® 0-2) + (a - /3) 0-3 (g) 0-3) , (75) 
where we have used the well-known Pauli matrix decomposition of the Bell states (see, e.g.. 



Ref. [10|). 



In order to determine the HS measure for the entangled two-parameter states p^"]^ we 
have to find the nearest separable states, which is usually the most difficult task to perform 
in this context. In Ref. [ll| a lemma is presented to check if a particular separable state is 
indeed the nearest separable state to a given entangled one: 
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\{\iij-){iij-\ + \iij-){r\) 



FIG. 2: Illustration of the two-qubit states pa 



1). 



Lemma 1. A state p is equal to the nearest separable state po if (^iT'd only if the operator 



C 



P - Pent - (P, P - Pent) 1 
Hp - Pent II 



(76) 



is an entanglement witness. 



Lemma [T] is used here in the following way. First, we calculate the separable state that 
has the nearest Euclidean distance in the geometric picture (Fig. [2]) and call this state p. 
But since the regarded picture does not represent the full state space (e.g., states containing 
terms like Oj a* (S> 1 or foj 1 ® cr* are not contained on the picture), we have to use Lemma [H 
to check if the estimated state p is indeed the nearest separable state po. 



1. Region I 

Let us consider first the entangled states located in the triangle region that includes the 
Bell state 10"*") , i.e. Region I in Fig. [2l For P = the states represent the isotropic state 
( !39|) . thus all results have to agree in this case with Eqs. (l63|) - (!66l) for d = 2. An entangled 
state in Region I is characterized by points, i.e. by the parameter pair {a,j3), constrained by 

a < p + a < + a > + (77) 
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The point in the separable region of Fig. [2] that is nearest (in the Euchdean sense) to the 
point is given by (| + which corresponds to the state 

^ f 1 + . ,1-2(3 \ 
Pp = 4 I 1 ^ i^i (S) ai - a2 (S) as) H — 0-3 ® 0-3 J . (78) 

For the difference of nearest-separable and entangled state we obtain 

h-pT!, = (79) 

where S is defined by 

S := cTi ® (Ti - 0-2 ® 0-2 + 0-3 (g) 0-3 . (80) 
Using the norm = 2^/3 we gain the HS distance 

\\p,-pT^\\ = ^{^~l-lf^)- (81) 

To check whether the state coincides with the nearest separable state po;/3 in fhe sense 
of the HS measure of entanglement (156|) (which has to take into account the whole set of 
separable states), we have to test — according to Lemma [1] — whether the operator 

is an entanglement witness. Remember that any entanglement witness A that detects the 
entanglement of a state p^'^* has to satisfy the inequalities ([5] 
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We calculate 



(P/3, PP - & = Tr pXp, - p-^) = (« - ^ - 1^) (83) 
and use Eqs. ( 1791) and (IHT]) to determine the operator C for the considered case, 

C^^^il-^). (84) 

Then we find 

(pS,^ = < 0, (85) 



3 3' ^ 

since the entangled states in the considered Region I satisfy the constraint a > |/3 + | . 
Thus the first condition of inequalities (IFTI) is fulfilled. 

Actually, condition fl85|) is just a consistency check for the correct calculation of operator 
C since by construction of C we always have (p™*,C) = — ||p — p'^'^^W < 0. Thus more 
important is the test of the second condition of inequalities fl57j) and in order to do it we 
need the following lemma: 

Lemma 2. For any Hermitian operator C that is of the form 

C = a (1 + Ci ((Ta; ® CTj; — (Tj, ® <Jy) + C2 O"^ <7z) CL ^ K-^, " 1 < Q < 1 (86) 

i/ie expectation value for all separable states is positive, 

(p, C) > Vp G 5 . (87) 

Proof. Any separable state p is a convex combination of product states and thus a 
separable two-qubit state can be written as the Bloch vector (see Refs. fiol. [lll|) 

with nf,m^' e IR, jn'"! < 1 , |m''| < 1 , Pfc > 0, ^^.p^ = 1 . (88) 
Performing the trace we obtain 

(p,C) = TipC = Y,Pka{l + ci [nlml - nlml) + C2 n>3 . (89) 

k 

We have 



ci (n^m^ - njmj) +C2n>^| < |n^| |m^| + |nj| |mj| + |m^| < 1, (90) 

{p,C) > Vp G ^. □ (91) 



y y 

and therefore 



Since the operator C f lM|) is of the form fISB]) we can use Lemma [2] to verify 

{p,C) > Vp G ^. (92) 
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Therefore C (18^ is indeed an entanglement witness and pp is the nearest separable state 
Pl3 = Po;/3 for the entangled states p™^ in Region I. 



Finally, we find for the HS measure of the states in Region I 

73/ 11 

— \ a 

2 1 3 3 



D{pt,) = IIPo;/3 - pZW = ^ ( « - i - i/^ ) • (93) 



2. Region II 

It remains to determine the HS measure for the entangled states p^"^ located in the 
triangle region that includes the Bell state \4>~) , i.e. Region II in Fig. [2l Here the entangled 
states are characterized by points {a, j3) , where the parameters are constrained by 

a < p + a > \(3-\, a < -f3-l. (94) 

The states in the separable region of Fig. [2] that are nearest to the entangled states 
in Region II are called pa,i3 and characterized by the points 



a\ _ f 1/3 (-l + 2a-/3) 
P J ~ I 1/3 (-2-2a + /5) 



(95) 



The necessary quantities for calculating the operator C are the following 

Pa,f3 - pTj3 = {a + l + p) {a^ ® cr^ - (g) cTj^ - 0-^ (g) (x^) , (96) 

l|Pa,/3-p5ll = ^(-«-l-/3) , (97) 

{Pa,f3, Ptp - Pa,^) = ^ (« + 1 + /?) , (98) 

so that C is expressed by 

C = (1 + ai ® ai - (T2 ® as - as ® (T3) . (99) 

2v 3 

To test C for being an entanglement witness we need to check the first condition of inequal- 
ities ( 1571) : we get 

(^-^,0 = -1= (« + ! + /?)< (100) 

as expected. Since operator C f l99|) is of the form fl86|) we apply Lemma [2] and obtain for 
the separable states 

(p,C') > VpG5. (101) 

Therefore also in Region II operator C f l99l) is indeed an entanglement witness and Pa,i3 is 
the nearest separable state p^,/? = Po;a/3 for the entangled states p'^^^^. 

For the HS measure of the states in Region II we find 
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C. Two parameter entangled states — qutrits 

The procedure of determining the geometry of separable and entangled states discussed 
in Sec. IVIIBI can be generalized to higher dimensions, e.g. for two-qutrit states. Let us 
first notice how to generalize the concept of a maximally entangled Bell basis to higher 
dimensions. A basis of maximally entangled two-qudit states can be attained by starting 
with a maximally entangled qudit state |0o) and constructing the other d? — 1 states in the 
following way: 

10,) = Ui®l\<po) ^ = l,2,...,c?2-l, (103) 

where {Ui} represents an orthogonal matrix basis ([1]) of unitary matrices and Uq usually 
denotes the unity matrix 1 (see Refs. 18, 19|). 



A reasonable choice to start with is the maximally entangled state fHOj) and using the 
WOE (see Sec. |V]) which is an orthogonal basis of unitary matrices. Such a construction has 



been proposed in Ref. [I6j. Then we set up the following d projectors onto the maximally 



entangled states - the Bell states: 

Pnk ■■= (f/nfc®l)lO«|(f/I,®l) n,k = 0,l,...,d-l. (104) 

We can express the Bell projectors as Bloch vectors by using the Bloch vector form ( l53l) of 
-Poo '■= 10+) (0+1 and the relations (indices have to be taken modd) 

Ulm = e^""^f/_„_,, (105) 

UnmUlk = e^™'f/„+i„+fc. (106) 



It provides for the Bell projector the Bloch form 



^ d-l 

= ;^Y1 Ulm ® U.l^ . (107) 



rf2 

m.,l=0 



In case of qutrits {d = 3) the 9 Bell projectors (11071) form an 8-dimensional simplex 
which is the higher dimensional analogue of a 3-dimensional simplex, the tetrahedron for 
qubits, see Fig. [31 This 8-dimensional simplex has a very interesting geometry concerning 
separability and entanglement (see Refs. 13, 30|). Due to its high symmetry inside — named 



therefore the magic simplex by the authors of Ref. [17| — it is enough to consider certain 
mixtures of Bell states which form equivalent classes concerning their geometry. 

We are interested in the following two-parameter states of two-qutrits as a generalization 
of the qubit case, Eq. (172|) . 



Pa,p = ^ g ^ 1 + aPoo + 13^ (Pio + P20) . (108) 



According to Ref. [17[ the Bell states represent points in a discrete phase space. The indices 
n, k of the Bell states can be interpreted as "quantized" position coordinate and momentum, 
respectively. The Bell states Pqq, Piq and P20 lie on a line in this phase space picture of the 
maximally entangled states, they exhibit the same geometry as other lines since each line 
can be transformed into another one. 
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FIG. 4: Illustration of the two-qutrit states pa^p (jlOSp . 



Inserting the Bloch vector form of PqO) -Pio and P20 (HO 71) we find the Bloch vector expan- 
sion of the two-parameter states (11081) 

Pa,/3 = ^fl + f«-fVi + (« + /5)f^2l , (109) 



9 V V 2, 
where we defined 

Ui := f/01 ® t/01 + f^02 ® t^02 + f/ii ® f/-ii + ® f/-i2 + ® f/-2i + f/22 ® , 
t/2 := f/10 ® U-io + f/20 ® f/-20 • (110) 

The constraints for the positivity requirement (po,/3 > 0) are 

«<^/?+l, «<-/? + !, «>^_1, (111) 

and for the PPT 

a < -p--, a > -p--, a < ^ + -. (112) 
- p 2' - 4^ 2' - 8 4 ^ ^ 

The Euclidean picture representing the HS space geometry of states (I108P is shown in Fig. IH 
where the parameter coordinate axes are non-orthogonal since in HS space they have to be 
orthogonal to the boundary lines of the positivity region, « = f — | and a = |/? + 1. It is 
shown in Ref. T3| that the PPT states pa,(3 are all separable states, so there are no bound 
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entangled PPT states of the form (11081) . To find the HS measure for the entangled two- 
parameter 2-qutrit states we apply the same procedure as in Sec. IVIIBI We determine the 
states that are the nearest separable ones in the Euclidean sense of Fig. [Hand use Lemma [1] 
to check whether these are indeed the nearest separable ones with respect to the whole state 



space (for other approaches see, e.g., Refs. |3ll. l32[|). 



1. Region I 

First we consider Region I in Fig. HI i.e., the triangle region of entangled states around 
the a-axis, constrained by the parameter values 

a < ^/5 + l, a<-/? + l, «>^ + ^- (113) 



In the Euclidean picture the point that is nearest to point {ot-,13) in this region is given by 
+ 



^ + ^/3,/3), which corresponds to the separable two-qutrit state 



with Ui and U2 defined in Eq. f lTTOD . 

For the difference of nearest-separable and entangled state we find 

where U = Ui + U2 (and is defined in Eq. (15^ ). and using for the norm \\U\\ = rfVs = 6^2 
we gain the HS distance 

IIp/^-pSII = • (116) 

It remains to calculate 

{Pl3, P/3 - Pa,p) = Tr p^(p^ - pa,p) = " ^ + (H^) 

to set up the operator 

\\pp-pT,f3\\ 6V2 

We test now whether it represents an entanglement witness, i.e., whether C (IllSp satisfies 
the inequahties ( 1571) . As expected we find 



To check the second condition of inequalities (1571) we set up the following lemma, similar to 
Lemma [2) 
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Lemma 3. For any Hermitian operator C that is of the form 

C = a{21 + CiUi + C2U2) a G R+, -1 < Q < 1 (120) 
the expectation value for all separable states is positive, 

(p,C) > \/peS. (121) 
Proof. Any separable two-qutrit state can be written as the Bloch vector 10, 111 

d-l d-l 



n,m=0 



Lk=0 



d-l 



n,m,l,k=0 

with n^^, m'l^eC, \n^\<l,\m^\<l, Pfc > 0, J^p^ = 1 



(122) 



where we define \n'^\'^ := Yjnm.Kim^r, 
Performing the trace we obtain 



(p,C) = Trptc = Y,Pk I 2a I 1 + ciJ^n.V + J] 



1=1 



with 



(123) 



{ni,n2,n3,ni,n5,ne,nT,n8) = (nQ^, ^02, <i, ^12, ^21. ^22. '^lO' ^^20) and 
(mi, m2, ms, m4, mg, me, ruj, mg) = (m*i, m*2, ml^i, mli2' "^-21^ "^-22^ "^-ic "^-20) • 

(124) 



We have 



and therefore 



i=l 



J=7 



(125) 
(126) 



{p,C) > Vp G 5. □ 
Since the operator C (11181) is of the form (11201) we can use Lemma [3] to verify 

(p,C') > ypeS. (127) 

Thus C (11181) is indeed an entanglement witness and p^ is the nearest separable state 
P/3 = Po;/3 for the entangled states p^"^^ in Region L 

For the HS measure of the entangled two-parameter two-qutrit states (11081) we find 



3 V 4 8'^' 



;i28) 
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2. Region II 

In Region II of Fig. |4] the entangled two-parameter two-qutrit states are constrained by 
a<^/3-^, a>^/5-^, a < -/3 + 1 . (129) 

The points that have minimal Euclidean distance to the points (a,/?) located in this region 
are characterized by 

(130) 



( ^ \ = ( (-2 + 20a + 5/3) \ 
[pj [ l/6(2 + 4a + /?) J 



and correspond to the states pa,i3- The quantities needed for calculating C are 

P«,/5-p5 = -;^(4a + 2-5/?)([/i-f/2) , (131) 

IIPa,/3-p5ll = ^(-4a-2 + 5/?) , (132) 

{Pa,p,Pa%-Pa,p) = ^ (4a + 2 - 5/?) , (133) 
so that operator C is expressed by 

C = ^ (21 + f/i - U2) . (134) 
The check of the first condition (!57|) for an entanglement witness gives, unsurprisingly, 

(p-^, (7) = -i= (4a + 2 - 5/3) < , (135) 

since 4a < 5/3 — 2 , Eq. (11291) . For the second test we use the fact that operator C ( 11341) is 
of the form (11201) and thus, according to Lemma [3l we obtain 

(p, C") > Vp G 5 . (136) 

Therefore C (11341) is indeed an entanglement witness and the states Pa,i3 are the nearest 
separable ones pa,/3 = Po- to the entangled two-parameter states (11081) of Region II. 
Finally, for the HS measure of these states we obtain 

D{pT^) = ||po;«,/^ - pT}^\\ = ^ (-4a -2 + 5/3) . (137) 



VIII. CONCLUSION 



In this paper we present three different matrix bases which are quite useful to decompose 
density matrices for higher dimensional qudits. These are the generalized Gell-Mann matrix 
basis GGB, the polarization operator basis FOB, and the Weyl operator basis WOE. Each 
decomposition we identify with a vector, the so-called Bloch vector. 
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Considering just one-particle states we observe the following features: 
The GGB is easy to construct, the matrices correspond to the standard SU(N) generators 
{N = d), but in general (in d dimensions) it is rather unpractical to work with the diagonal 
GGM (E]) due to its more complicated definition. On the other hand, the Bloch vector itself 
has real components, which is advantageous, they can be expressed as expectation values of 
measurable quantities. For example, in 3 dimensions the Gell-Mann matrices are Hermitian 
and the Bloch vector components can be expressed by expectation values of spin 1 operators. 
The FOB is also easy to set up, all you need to know are the Clebsch-Gordan coefficients 
which you find tabulated in the literature. However, the Bloch vector contains complex 
components. The Weyl operators for the WOB are also simple to construct, they are non- 
Hermitian but unitary operators. The Bloch vector itself turns out very simple, however, 
with complex components. Let us note that in 2 dimensions all bases are equivalent since 
they correspond to Fauli matrices or linear combinations thereof. 

In case of two-qudits we have studied the isotropic states explicitly and find the following 
characteristics: 

In the GGB the Bloch vector fH7|) with expression fHSl) is more complicated to construct, 
in particular the diagonal part B (l45l) (see Appendix I A 21) . In the FOB the Bloch vector (152|) 
with expression (15T!) can be easily set up by the knowledge of the Clebsch-Gordon coefficient 
sum rule (1491) and in the WOB the Bloch vector (l55l) with definition is actually most 
easily to construct. 

The Hilbert-Schmidt measure of entanglement has been calculated explicitly for the 
isotropic two-qudit states and we want to emphasize its connection to the maximal violation 
of a generalized Bell inequality (Theorem [T]), an inequality for the entanglement witness. 

We demonstrate the geometry of separability and entanglement in case of qubits by 
choosing so-called two-parameter states, Eq. fl72l) . i.e., planes in the tetrahedron formed 
by the Bell states (see Fig. [3]). These states refiect already the underlying geometry of 
the Hilbert Space and they are chosen with regard to the description of qutrit states, a 
generalization into higher dimensions. To a given entangled state we determine the nearest 
separable state, calculate the corresponding entanglement witness and the Hilbert-Schmidt 
measure in the relevant Regions I and II (see Fig. [2]). 

In case of qutrits it is quite illustrative to demonstrate the geometry of separability and 
entanglement in terms of two-parameter states fllOSp . These states set up a plane in the 8- 
dimensional simplex formed by the 9 Bell states and are easy to construct within the WOB. 
Due to the high symmetry of the magic simplex we may restrict ourselves to a certain 
mixture of Bell states, Eq. (llOSp . which exhibits the same geometry as other lines. Within 
the WOB it is quite easy to find the Bloch vector form (I109P of the two-parameter states. 
It is straightforward to find for a given entangled state in the relevant Regions I and II (see 
Fig. H]) the nearest separable state and the corresponding entanglement witness. The easy 
calculation of the Hilbert-Schmidt measure of entanglement is a great advantage in this 
case and its result of high interest since it is quite difficult to calculate other entanglement 
measures for higher dimensional states, like the entanglement of formation. 

It turns out that the Weyl operator basis is optimal for all our calculations, the reason is 
that entanglement — the maximally entangled Bell states — is in fact easily constructed by 
unitary operators a la Weyl. 
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APPENDIX A 

1. Proof of Orthogonality of GGB 

We want to proof condition for the GGB which consists of the d"^ - 1 GGM (El), (HD, 
([5]) and the dx d unity 1. Since all GGM are Hermitian (thus TrA^Aj = TiA^Aj = TrAjAi) 
it suffices to proof the following conditions: 



TrAfA^'^ = 2(5^''"(5'=" (Al) 

TrAfA;^" = 2(5-'™5'=" (A2) 

TrA'A™ = 2 5'™ (A3) 

TrAfA^*^ = (A4) 

TrAfA'" = (A5) 

TrAfA'" = 0. (A6) 



Proof of condition (lAip . Inserting definition we have 

d 

TrAfAr = J2(^\{\j){k\ + \k){j\){\m){n\ + \n){m\)\l) 
1=1 

= 5^((/|j)(A;|m)(n|0 + {k\n) {m\l) + m{j\m){n\l) + {l\k){j\n){m\l)) 
I 

^jn^km _j_ ^jm^kn _j_ ^kn^jm _|_ ^km^jn 

= 2(5^'^5'=", (AT) 
where we used in the last step that S^"-S^"^ = since we have j < k and m < n. 

Proof of condition (lA2p . This case is equivalent to the one before apart from changed 
signs that do not matter 

TrA-'^A"*" = — -|- §3^^^^ _|_ §kn^jm _ ^km^jn 

= 2 5^'"5'=". (A8) 
Proof of condition flA3l) . Using definition ([5]) and denoting 



(A9) 
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where / < m without loss of generahty, we get 

d I m 

TrA'A'" = CiCra^(^^^{p\k){k\n){n\p) + lm{p\l + l){l + l\m + l){m + l\p) 

p=l k=l n=l 

m ' 

k=l n=l 

m 



.^{p\k){k\m + l){m + l\p) - l^{p\l + I) {I + l\n) {n\p) 

k=l n=l 

/ I m \ 

CiCm [l + lm5''^ - m^(5^(™+^) - . (AlO) 



k=l n=l 

Using the fact that (J'^^'^+i) = for m > A; and 



= <'° '^^ = '^ (All) 

^ I / if / < m ^ ^ 



n=l 

we obtain 



TrA'A"^ = + 1)5'" = 2 5'"^. (A12) 

Proof of condition ( ]A4I) . Analogously to the proofs (1A7|) and (1A8|) we find 

Tr Af A™" = z (-^J'^^^™ + - ^J'"^*^" + ^i"^^™) = o. (A13) 

Proof of condition f lASp . Inserting definitions ([3]) and (jSj) gives 

TrA^A™ = C^y^ f - m(p|A;)(j|m + l)(m + l\p) - ni{p\j){k\m + l)(m + l|j9) 
p=i 

m m 
n=l n=l 



= - 2m + 2 ^ ^'^'(J^' 

1=1 

= 0, (AM) 
since per definition we have j < k . 



Proof of condition flA6p . This proof is equivalent to the previous one since constant 
factors in front of the terms do not matter. 



2. Calculation of term B in GGB 



To obtain the Bloch vector notation of term B fH51) we insert the standard matrix ex- 
pansion (HM for the case j = k. We split the tensor products in the following way 

B = i ( 5i + ^2 + ^3 + ^4 + i 1 ® 1 j , (A15) 
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where the terms _Bi, . . . , S4 are introduced by (note that A° = 0) 

. T d-j-1 
j=l V n(=0=0 



5i = V — -A^~' ® + V -A^+" ® A^+" (A16) 

lf\ 2j A. 2(j+n)(i + n + l) ^ ^ 



^2 = y I -y\l—, — — -A^'-i®A^-+' 



d-j-i 



J -I 



^ y 4j(j+n)(j+n + l) 

+ V — ® A^'+' (A17) 

n^W=o V(J +n){j+n + + + / + 1) ' 

^3 = ly" ( -j:^A^-i® 1 + y , ^ A^'+" 1 ) (MB) 



^4 = -,y(-\- — -l^A'-' + V , ^ 1 ^ A^'+' 1 . (A19) 



d 

Only the first term Bi (IAT6D Gfives a contribution 



Si = ( ^ ^ + ^ ) A™ ® A™ = - V A"^ ® , (A20) 

^ \2(m + 1) 2m(m + 1) J 2^ ' ^ ^ 

m=l ^ ^ ' ^ ' ' m=l 



^ \2(m + 1) 2m(m + 1) 
whereas the remaining terms vanish: 

d-l 



m<p, m,p=l 



B2 = y -J-, — — , + , ^ A-^A^ 



m>p, m,p=l 



+ y I-a/tT TT^T TT+ , ^ )a-0A^ 



4(p + l)m(m+l) ^4p(p + l)m(m + 1) 



> — — + y — — A™ (g) A*' 

^m<p 2Vm(m + l)p{p + 1) 2v/m(m + l)p{p + 1) J 

= 0, (A21) 
and in quite the same manner 

, d-l 

53 = - V , A™ ® 1 = , 

= -y , 1 ® AP = . (A22) 



Thus we find the following Bloch vector of B (1431) 



1 "^"^ 1 
B = _^A™®A™ + -1®1. (A23) 

m=l 
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3. Proof of Orthonormality of WOB 

For proofs relevant in the WOB we often need the equivalence 

n=0 k 

So we use Eq. (lA24p to proof the orthonormality of the Weyl operators (130]) 

TrUl^Uij = ^ ^ e^^^'-'^") (p|(A; + m)modc/)(A;|fc)((fc + j)modc/|p) 

P=0 k,k=0 
d-1 d-1 

= ^ J] {p\ {k + m) mod rf) ((A; + j) mod 6^-^ 

P=0 k,k=0 
k=0 



= d5ni5mj- (A25) 

4. Expansion into WOB 

Formula fl53l) for the Bell state in terms of WOB we derive in the following way. We 
express the standard matrices by the WOB fl36l) . rewrite the indices and separate the non- 
vanishing terms 

j,k=i 

d-1 d-1 

^ ^ ^ Ul(k-j)modd ® Uv{k-j)modd 

j,k=Ol,l'=0 



Ci3 

m,fc=0 /,/'=0 



V m m l,l'-l+l'=d 

+ i E E f E ^'"^ u,^ ® f/,^ . (A26) 



The last term in Eq. (1A26|) vanishes due to relation (lA24p . Identifying Uqo = 1 and using 



the notation with negative values of the index /, which have to be considered as mod d , we 
gain the formula 

1 1 '^'^ 

k+>«l = ^1®! + ^ E^'™®^"'™' (/,"^) 7^ (0,0). (A27) 

Z,m=0 
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